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Abstract 

The integrated Sachs-Wolfe (ISW) effect in a A dominated uni¬ 
verse can be an important factor in the evolution of cosmic microwave 
background fluctuations. With the inclusion of cosmological constant 
we present the complete analytic solution of the covariant linear per¬ 
turbations of the Einstein equations in the Newtonian gauge, in the 
case of a spatially flat (K = 0), Friedman-Robertson-Walker (FRW) 
universe, filled with pressureless ideal fluid. We use the analytic time 
dependence of the perturbation potentials to derive the anisotropy 
power spectrum of the late ISW effect. We choose the scale invariant 
Harrison-Zeldovich spectrum for obtaining the Cg momenta and show 
the agreement of our results with earlier, numeric calculations. 
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1 Introduction 

It is well known that several effects give contribution to the fluctuations of 
the cosmic microwave background radiation (CMBR). It is usual to divide 
these contributions to primary and secondary anisotropies. The dominant 
part of primary anisotropies comes from the Sachs-Wolfe effect I]. In a A 
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or a curvature dominated universe the gravitational potentials of the pertur¬ 
bations are varying in time, giving rise to the appearance of the ISW effect. 
In linear theory for Oo = 1, flat universe, the potentials are constants so this 
appearance does not occur. 

The ISW effect in a A dominated universe has been studied first by [2] 
and recently in more detailes, as different measurements e.g. Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAP) [2], la type supernovae observations J3] 
and the Sloan Digital Sky Survey (SDSS) s.5j imply the accelerating expansion 
of the universe. In the light of this unexpected result the understanding of 
the properties of A models becomes more and more important. In a thorough 
overview [ 6 | Hu and Sugiyama summarize the CMBR anisotropies concern¬ 
ing many different scenarios. In |7 Hu and White investigate the CMBR 
anisotropies in the weakly coupled regime (i.e. when the photons travel 
freely from the last scattering surface to present) including the A dominated 
epoch. These works give an accurate picture of present understanding of 
the microwave background fluctuations using both analytic and numerical 
methods, although none of them have the general analytic expression for the 
perturbation potentials in the presence of cosmological constant. 

In our paper, following the treatment of earlier works we calculate the late 
ISW effect in a A dominated universe. In Sec. 2, we present the metric with 
linear perturbations for a flat, homogeneous and isotropic universe. Using 
the analytic results of [ 8 ] we obtain the complete solution of the perturbed 
problem in the Newtonian gauge for scalar perturbations. In Sec. 3, we 
calculate the Ci momenta and the power spectrum of the anisotropies for 
the late ISW effect and show that our results are in good agreement with 
earlier full numeric calculations. 


2 The Newtonian potentials of the perturba¬ 
tions 

The unperturbed metric in a homogeneous and isotropic universe is the FRW 
metric 


ds 2 = — dt 2 + a 2 (t)^ llu dx fJ 'dx l/ , ( 1 ) 

with the three metric 7 ^ of a space with constant spatial curvature K. The 
scale factor a{t ) can be expressed in the case of a flat (K = 0) space and 
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pressureless (p = 0) matter source as [.9] 


a = a 0 sinh 2//3 (C't + C 0 ) , pa 3 = Cm , a 0 = 



where, and throughout this paper we use units in which the gravitational 
constant G = 1/87r, the speed of light c = 1 and we set C 0 = 0. 

The linear perturbations of the metric tensor for a spatially flat, FRW 
universe have the general form 



( 3 ) 


9a/3 — a? [(1 + 2 H l ) r ) a 0 + 2 H Ta/3 \ , 


where the functions A, Hi, B a and Hi a p give a complete representation 
of the metric and H^ap is a 3 x 3 trace-free tensor. In Bardeen’s covariant 
linear perturbation formalism the Einstein equations can be decoupled into a 
set of ordinary differential equations by employing scalar, vector and tensor 
eigenmodes of the Laplacian operator, which form a complete set. In the 
following we will represent the perturbation quantities with the amplitudes 
of the eigenmodes corresponding to the fcth wavenumber, see QDj. 

In covariant linear formalism three independent, gauge invariant quanti¬ 
ties can be constructed, namely T and $ from the metric tensor amplitudes 
and V from the matter velocity ma. These gauge invariant quantities are 



( 4 ) 


where a dot denotes the derivative with respect to the comoving time coor¬ 
dinate t and v is the velocity perturbation of the matter held. We use the 
complete first order solution given in comoving coordinates [S|, to obtain the 
metric tensor amplitudes in the Newtonian (or longitudinal) gauge. This 
transformation can be done with the help of the gauge invariant quantities 
in Eq. 0 using the Newtonian gauge conditions H^ = B N = 0, where 
the superscript N refers the Newtonian gauge. Performing the calculations 
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yields the solution 


T = A N 

V = v N 


C cosh(Ct) {3 C M B 0 (k) - 2Ca§ [ATi(Jfc) 
3 a 0 k 2 CM sirili 5 ' 73 (C'^) 

-v& , 


K2(k)I(t)]} t (5) 

( 6 ) 


1 f Co 4 
sinh 1 ^ 3 (Ct) { kC\f 
2V 3 Cag sinh 1/3 (C't) 
kCj^ cosh(Ct) 


[Ad(A:)-/l 2 (fc)/(t)] 
-K 2 (k) . 


3B 0 (k) \ 
2 kal J 


(7) 


The density perturbation from the held equations also can be obtained im¬ 
mediately as 


Sp N = 


cosh (Ct) ( 4agC 2 
sinh 3 (Ct) [ 3 C M 


[Ki(k) - K 2 {k)I{t)) 


^fB 0 (k) 


a; 


o 


+ 


3 sinh 4/3 (Ct) 
k 2 cosh(Ct) 


u2 o4/3p2 3 

^H 0 (k) + ' ° lf 2 (fc) 


c 


M 


cto sinh 7 ^ 3 (C't) _3 Cm Cap 
Using the gauge restrictions of jHj 


5 0 (fc) . 


Ki(k) 

H 0 (k) 


k 2 C 


M 


2 On 


At 



tf 2 (fc) 

5 0 (fc) 


3fc 2 UM 
2 7 ^alC 2 k ’ 


0 , 


( 8 ) 

(9) 


we can eliminate the non-physical (gauge) terms from the solutions. It is 
shown in m that the Ki(k) term in Eq. (JHJ) corresponds to relatively de¬ 
creasing density perturbations. Calculating the anisotropy power spectrum 
of the late ISW effect we are interested in the relatively growing modes, thus 
in the following we keep only the K 2 (k ) term in the T and <f> potentials. 
Applying these restrictions one can rewrite the solution as 


<f> 


_ 1 cosh (Ct) 

2 4 / 3 sinh 5 / s (Ct) W *’ 
-T , 
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( 10 ) 

( 11 ) 





















V 


( 12 ) 


3k [ I(t) _ sinh 1/3 (Ct) 
4 Ca 0 _2 1 / 3 smh 4 / 3 (Ct) cosh (Ct) 

k 2 cosh (Ct)I(t) 

P ~ 2 1 /3flg sinh 3 (Ct) k ’ 


(13) 


where 


I(x) 

x{t) 


2^6^ +V3i ^E(x) 

^ = sinh 2 / 3 (Ct) , 
a 0 


3 + y/3 i 
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nl/s V^i 1 ~ X ) 

\/x 3 + 1 


,( 14 ) 


(15) 


and E, F are the following incomplete elliptic integrals sm 


E(x) 


F(x) 


/ 

(3 + V3i)x 

i-V3\ 

h 

2(i + l) 

2 J 

/ 

(3 + a/3 i)x 


U 

2(x + l) 

2 J 


(16) 


Throughout this section we have used the comoving time coordinate t 
in our calculations. In order to derive the anisotropy power spectrum it 
is advisable to introduce the conformal time variable with the definition 
adrj = dt. To transform our results we replace x(t) with x(rj), where 


x{rj) 


-1 + 


3 + y/3i 


3 + V3i-2sn 2 U 1 ^ m - ) Ca 0 r) 


i-V3 

2 


(17) 


and sn (rj,z) is the Jacobi elliptic function (cf. Appendix). 


3 The power spectrum of the late ISW effect 

In the previous section we obtained the complete analytic solution of the 
scalar perturbation quantities in the Newtonian gauge. In possesion of this 
result we can deduce the anisotropy power spectrum of the late ISW effect. 
The method of the calculation is well known jHJ, therefore we do not present 
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it in details, and hereafter we discuss only the main steps and the most 
important equations of the procedure. 

The power in the £th multipole is denoted by Cg, where 

21+l c = V_ [°° dk 3 \e e (r},k)\ 2 
47t e 27r 2 Jo k 2^ + 1 



and @£ is the multipole decomposition of the fractional temperature fluctua¬ 
tion |.6j. In the A dominated universe the perturbation potentials vary with 
time and photons experience differential redshifts due to the gradient of T, 
which do not yield equal and opposite contributions as the photons enter and 
exit the potential well (as is the case in the ordinary Sachs-Wolfe effect), and 
time dilation from <f>. The sum of these contributions along the line of sight 
is called the ISW effect. The total contribution is 


= [ 0 o + k)ji[k(rjo - 77 *)] (19) 

rvo 

+ / [+' - $'] ( 77 , k)j e [k(rj 0 - r])]dr] , 

Jr]* 

where the first term represents the ordinary, while the second term the in¬ 
tegrated Sachs-Wolfe effect. In Eq. (fill je is the Bessel function, 77 * and 770 
are the values of the conformal time at decoupling and present respectively, 
and a prime denotes the derivative with respect to 77 . I 11 the following we 
will focus on the contribution of the ISW term to the anisotropy power spec¬ 
trum of CMBR. It is usual to divide the ISW effect into early and late parts. 
The early ISW effect corresponds to the variations in the perturbation po¬ 
tentials which occurs in the time interval between matter-radiation equality 
and recombination, while the late ISW effect corresponds to the variations 
caused by the presence of the cosmological constant. In the A dominated 
epoch (which occurs recently) for adiabatic perturbations this late ISW term 
in Eq. PI) can be approximated with the formula Q 

P10 Tn 

j [* - *']jilKm - v)]dv - ■fl'S' - *'],=„ . (20) 

where rjk = 770 — (£ + 1 / 2 )//c and the A integral is given by 


h 



V? r|2(7 + i)] 

2 r[i (7 + 2)] ' 


( 21 ) 
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Substituting the results above into Eq. m the C( momenta of the late ISW 
effect have the form 


r dsw 


2V 


/ r[(£+1)/2] \ 
\r[(< + 2)/2]/ 


roo 

/ 'b /2 (%, k)dk . 
Jo 


( 22 ) 


In A cosmological models the perturbation potential T is related to the 
growth factor D as T ~ D/a, where 


D(a) ~ H 


da 

(Ha) 3 


(23) 


and H is the Hubble parameter dehned as H = a'/a 2 . The integral in Eq. 
© has been evaluated so far only with numerical treatment in the presence 
of cosmological constant |> • With the solution of the perturbation potential 
\b, we can express the growth factor analytically as follows 


D ~ aT ~ coth (Ct)I(t) 


irr) 



(24) 


Inserting the solution of T into Eq. (12211 the argument of the integral becomes 

2 

l*(0, k)\\ (25) 

where \b(0, k) = —314/10. 

In order to calculate the C/ SH momenta, we choose a general power law 
initial spectrum V(k) = fc 3 |'b(0, k)\ 2 = Bk ^ 1 , where B is a constant. It is 
usual to assume the scale invariant Harrison-Zeldovich spectrum with n — 1. 
The matter and dark energy densities of the universe are set to be Ho = 0.3 
and Ha = 0.7 respectively, according to the combined results of the WMAP 
and SDSS surveys. For the Hubble parameter H 0 — h ■ 100 km/s/Mpc, we 
use h — 0.7. 

Having in hand the solution of the perturbation potentials we have ex¬ 
pressed the argument of the integral in Eq. (USD analytically, although the 
integration over k has still to be done numerically. After performing the 
calculations the obtained power spectrum of the late ISW effect is plotted 
on Fig. 1. These results are in good agreement with earlier calculations, 
where the explicit analytic form of the \b potential has not been known in 
the presence of cosmological constant mm®. 


vp'2 = ^2 ~ 2/3 C 
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Jx( 1 + x 3 ) 3a;3 
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Figure 1: The late ISW effect in a A dominated universe 

4 Concluding remarks 

In the present paper we calculated the anisotropy power spectrum of the late 
ISW effect in a A dominated universe. We obtained the complete analytic 
solution for the growth factor and the scalar quantities of the perturbations 
using the Newtonian gauge in the presence of cosmological constant. We 
also expressed the argument of the ISW integral analytically for an arbitrary 
initial power spectrum of the fluctuations. We showed the agreement of our 
results with earlier numerical calculations in the case of the scale invariant 
Harrison-Zeldovich spectrum. 

Our results are relevant in calculating the anisotropy power spectrum in 
the presence of cosmological constant. In the case of an arbitrary initial spec¬ 
trum of the fluctuations the anisotropy spectra can be more easily calculated 
with the help of the obtained analytic formulae. 
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5 Appendix 


The time dependence of the scale factor on the comoving time t is given by 
Eq. ((21) and Eq. (fTT)|) as 


a = a 0 sinh 2,/3 (Ct) = a 0 x . 

Using the definition of the conformal time 

adi] = dt , 

7] can be expressed with the following integral of the variable x 


1 f dt 

7i = — — = 


dx 


a o J x 2Ca 0 J x{ 1 + x 3 ) 
Eq. (I2%1) can be integrated analytically with the solution 


(26) 


(27) 


(28) 


7](x) = 


j~3 \J3 — V3i F 

( 

(3 + \/3 i)x 

i-V 3 ] 

/ 2 Ca 0 


2(i+l) 

2 J 


(29) 


where F is the same incomplete elliptic integral given in (d. Taking its 
inverse, one can obtain the conformal time dependence of x as 


x (v) = -i + 


3 + \/3i 


3 + y/3i - 2 sn 2 U^j^Caori 


i-\J 3 

2 


(30) 


with the Jacobi elliptic function sn (r],z). In spite of the fact that rj(x) and 
x(t]) have complex arguments, both functions are real valued. 
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